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A  model  of  thermoelectric  generator-driven  thermoelectric  refrigerator  with  external  heat  transfer  is  pro¬ 
posed.  The  performance  of  the  combined  thermoelectric  refrigerator  device  obeying  Newton’s  heat 
transfer  law  is  analyzed  using  the  combination  of  finite  time  thermodynamics  and  non-equilibrium  ther¬ 
modynamics.  Two  analytical  formulae  for  cooling  load  vs.  working  electrical  current,  and  the  coefficient 
of  performance  (COP)  vs.  working  electrical  current,  are  derived.  For  a  fixed  total  heat  transfer  surface 
area  of  four  heat  exchangers,  the  allocations  of  the  heat  transfer  surface  area  among  the  four  heat 
exchangers  are  optimized  for  maximizing  the  cooling  load  and  the  coefficient  of  performance  (COP)  of 
the  combined  thermoelectric  refrigerator  device.  For  a  fixed  total  number  of  thermoelectric  elements, 
the  ratio  of  number  of  thermoelectric  elements  of  the  generator  to  the  total  number  of  thermoelectric 
elements  is  also  optimized  for  maximizing  both  the  cooling  load  and  the  COP  of  the  combined  thermo¬ 
electric  refrigerator  device.  The  influences  of  thermoelectric  element  allocation  and  heat  transfer  area 
allocation  are  analyzed  by  detailed  numerical  examples.  Optimum  working  electrical  current  for  maxi¬ 
mum  cooling  load  and  COP  at  different  total  number  of  thermoelectric  elements  and  different  total  heat 
transfer  area  are  obtained,  respectively. 
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1.  Introduction 

Semiconductor  thermoelectric  power  generation,  based  on  the 
Seebeck  effect,  and  semiconductor  thermoelectric  cooling,  based 
on  the  Peltier  effect,  have  very  interesting  capabilities  with 
respect  to  conventional  power  generation  and  cooling  systems 
[1-3].  The  absence  of  moving  components  results  in  an  increase 
of  reliability,  a  reduction  of  maintenance,  and  an  increase  of  sys¬ 
tem  life;  the  modularity  allows  for  application  in  a  wide-scale 
range  without  significant  losses  in  performance;  the  absence  of 
a  working  fluid  avoids  environmental  dangerous  leakage;  and 
the  noise  reduction  appears  also  to  be  an  important  feature. 
Thermoelectric  generator  and  refrigerator  have  been  used  in 
military,  aerospace,  instrument,  and  industrial  or  commercial 
products,  as  a  power  generation  or  cooling  devices  for  specific 
purposes.  Many  researchers  are  concerned  about  the  physical 
properties  of  thermoelectric  material  and  the  manufacturing 
technique  of  thermoelectric  modules.  In  addition  to  the  improve¬ 
ment  of  the  thermoelectric  material  and  module,  the  system 
analysis  and  optimization  of  thermoelectric  generator  and  refrig¬ 
erator  are  equally  important  in  designing  high-performance 
thermoelectric  generators  and  refrigerators. 
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In  general,  the  conventional  non-equilibrium  thermodynamics 
[1,4]  is  used  to  analyze  the  performance  of  single-stage  one-  or 
multiple-element  thermoelectric  generators  [5-13]  and  refrigera¬ 
tors  [14-25],  Non-equilibrium  thermodynamics  is  a  progress  of 
classical  thermodynamics,  it  considers  some  specific  phenomena 
such  as  Seebeck  effect  and  Peltier  effect.  All  objects  of  these 
researches  are  independent  thermoelectric  devices,  that  is,  they 
generate  direct-current  power  source  for  users  (thermoelectric 
generator)  or  need  a  direct-current  power  source  to  provide  direct 
current  (thermoelectric  refrigerator).  In  some  special  fields,  the 
heat  rejected  from  the  thermal  machine  may  drive  a  thermoelec¬ 
tric  refrigerator  through  the  use  of  a  thermoelectric  generator,  so 
that  the  thermoelectric  refrigerator  does  not  need  an  independent 
power  source.  Such  a  new  system  is  different  from  the  traditional 
thermoelectric  systems  merely  consisting  of  a  thermoelectric  gen¬ 
erator  and  a  refrigerator.  These  systems  dispense  with  complicated 
pipelines  and  heat  insulation,  so  they  can  be  used  in  many  special 
fields  such  as  aircraft  and  submarine.  Chen  et  al.  [26]  and  Khattab 
and  El  Shenawy  [27]  built  a  model  of  this  kind  of  combined  system, 
i.e.  single-stage  thermoelectric  refrigerator  driven  by  single-stage 
thermoelectric  generator,  and  analyzed  the  performance  of  the 
device. 

The  theory  of  finite  time  thermodynamics  or  entropy  generation 
minimization  [28-35]  is  a  powerful  tool  for  performance  analysis 
and  optimization  of  practice  thermodynamic  processes  and 
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Nomenclature 

A 

coefficient  of  system  stable  working  electrical  current 

Greek  symbols 

equation 

a 

Seebeck  coefficients  of  P-  and  N-  type  semiconductor 

F 

heat  transfer  surface  areas  of  heat  exchanger 

legs 

f 

heat  transfer  surface  area  ratio 

£ 

coefficient  of  performance  (COP)  of  combined  thermo¬ 

I 

K 

working  electrical  current  (A) 

thermal  conductance  of  a  semiconductor  couple  (W/K) 

electric  refrigerator  device 

k 

heat  transfer  coefficient  of  heat  exchanger  (W/K) 

Subscripts 

M 

total  number  of  thermoelectric  elements  pairs  of  whole 

1 

parameter  of  thermoelectric  generator 

device 

2 

parameter  of  thermoelectric  refrigerator 

m 

number  of  thermoelectric  elements  pairs  of  thermoelec¬ 

HI 

heat  source  of  thermoelectric  generator 

tric  generator 

H2 

heat  sink  of  thermoelectric  refrigerator 

n 

number  of  thermoelectric  elements  pairs  of  thermoelec¬ 

LI 

heat  sink  of  thermoelectric  generator 

tric  refrigerator 

L2 

heat  source  of  thermoelectric  refrigerator 

T 

temperature  (K) 

S 

practical  solution  of  working  electrical  current  equation 

Q. 

R 

rate  at  which  heat  is  transferred  (W) 

total  internal  electrical  resistance  of  a  semiconductor 

couple  (C2) 

ratio  of  number  of  thermoelectric  element  pairs  of  the 
thermoelectric  generator  to  total  number  of  thermoelec¬ 
tric  element  pairs  of  the  combined  irreversible  device 

opt 

optimum  parameter 

devices.  Finite  time  thermodynamics  or  entropy  generation  mini¬ 
mization  is  the  method  of  modeling  and  optimization  of  various 
thermodynamic  processes  and  devices  that  owe  their  thermody¬ 
namic  imperfection  to  heat  transfer,  mass  transfer,  and  fluid  flow 
and  other  transport  processes.  It  bridges  the  gaps  not  only  between 
thermodynamic  and  heat  transfer,  mass  transfer,  fluid  mechanics, 
and  other  transport  science,  but  also  between  the  physics  and  the 
engineering.  It  thermodynamically  optimize  performance  of  real  fi¬ 
nite-time  and/or  finite-size  thermodynamic  systems  with  the  irrev¬ 
ersibilities  of  heat  transfer,  fluid  flow  and  mass  transfer  toward 
decreasing  the  irreversibility  of  the  total  system.  Some  authors 
have  investigated  the  performance  of  thermoelectric  generators 
[36-48]  and  thermoelectric  refrigerator  [49-57]  using  the  combi¬ 
nation  of  finite  time  thermodynamics  and  non-equilibrium  ther¬ 
modynamics.  They  analyzed  the  effect  of  finite-rate  heat  transfer 
between  the  thermoelectric  device  and  its  external  heat  reservoirs 
on  the  performance  of  single-element  single-stage  thermoelectric 
generators  [36-43]  and  refrigerator  [49-52],  They  also  investigated 
the  characteristics  of  single-stage  multi-element  thermoelectric 
generators  [44-48]  and  thermoelectric  refrigerator  [53-57]  with 
the  irreversibility  of  finite  rate  heat  transfer,  Joulean  heat  inside 
the  thermoelectric  device,  and  the  heat  leak  through  the  thermo¬ 
electric  couple  leg.  However,  all  of  those  were  performed  only  for 
independent  thermoelectric  devices.  There  has  been  no  investiga¬ 
tion  concerning  the  performance  analysis  and  optimization  for  sin¬ 
gle-stage  thermoelectric  refrigerator  driven  by  single-stage 
thermoelectric  generator  published  in  the  open  literature. 

On  the  basis  of  the  exo-reversible  model  of  a  single-stage  ther¬ 
moelectric  refrigerator  driven  by  a  single-stage  thermoelectric 
generator  without  external  irreversibility  built  in  Refs.  [26,27],  a 
model  of  thermoelectric  generator-driven  thermoelectric  refrigera¬ 
tor  with  external  heat  transfer  is  built.  The  performance  of  the 
combined  thermoelectric  refrigerator  device  obeying  Newton’s 
heat  transfer  law  is  analyzed  using  the  combination  of  finite  time 
thermodynamics  and  non-equilibrium  thermodynamics.  Two  ana¬ 
lytical  formulae  for  cooling  load  vs.  working  electrical  current,  and 
the  coefficient  of  performance  (COP)  vs.  working  electrical  current, 
are  derived.  For  a  fixed  total  heat  transfer  surface  area  of  four  heat 
exchangers,  the  allocations  of  the  heat  transfer  surface  area  among 
the  four  heat  exchangers  are  optimized  for  maximizing  the  cooling 
load  and  the  coefficient  of  performance  (COP)  of  the  combined 
thermoelectric  refrigerator  device.  For  a  fixed  total  number  of 


thermoelectric  elements,  the  ratio  of  number  of  thermoelectric  ele¬ 
ments  of  the  generator  to  the  total  number  of  thermoelectric  ele¬ 
ments  is  also  optimized  for  maximizing  both  the  cooling  load 
and  the  COP  of  the  combined  thermoelectric  refrigerator  device. 
The  influences  of  thermoelectric  element  allocation  and  heat  trans¬ 
fer  area  allocation  are  analyzed  by  detailed  numerical  examples. 
Optimum  working  electrical  current  for  maximum  cooling  load 
and  COP  at  different  total  number  of  thermoelectric  elements 
and  different  total  heat  transfer  area  are  given,  respectively. 

2.  Model  of  a  thermoelectric  generator-driven  thermoelectric 
refrigerator  system 

A  schematic  diagram  of  a  thermoelectric  generator-driven 
thermoelectric  refrigerator  device  is  shown  in  Fig.  1.  The  device 
consists  of  an  irreversible  single-stage  multi-element  thermoelec¬ 
tric  generator  and  an  irreversible  single-stage  multi-element 
thermoelectric  refrigerator  in  series  with  internal  and  external 
irreversibilities.  The  direct-current  power  source  of  the  refrigerator 
is  the  direct-current  power  output  of  the  generator. 

The  irreversible  thermoelectric  generator  is  composed  of  m 
pairs  of  thermoelectric  elements.  Each  element  is  composed  of 
P-type  and  N-type  semiconductor  legs.  The  thermoelectric  power 
generation  element  is  assumed  to  be  insulated,  both  electrically 
and  thermally,  from  its  surroundings,  except  at  the  junction- 
reservoir  contacts.  The  internal  irreversibility  is  caused  by  Joulean 
electrical  resistive  loss  and  heat  conduction  loss  through  the  semi¬ 
conductor  between  the  hot  and  cold  junctions.  The  Joulean  loss 
generates  an  internal  heat  I2R,  where  R  is  the  total  internal  electri¬ 
cal  resistance  of  the  semiconductor  couple  and  /  is  the  electrical 
current  generating  from  the  semiconductor  couple.  The  conduction 
heat  loss  is  K(T'm  -  T[, ),  where  K  is  the  thermal  conductance  of  the 
semiconductor  couple,  T'm  is  the  hot  junction  temperature,  and  T'n 
is  the  cold  junction  temperature.  Finite  rate  heat  transfers,  i.e.  the 
temperature  differences  ( Tm  -  T'm)  and  (Tu  -  Tu),  where  Tm  and 
Tii  are  the  temperatures  of  the  heat  source  and  heat  sink  of  the 
thermoelectric  generator,  respectively,  cause  the  external  irrevers¬ 
ibility.  For  the  thermoelectric  generator,  the  rate  of  heat  transfer  at 
the  hot  junction  is  Qm,  and  the  rate  of  heat  transfer  at  the  cold 
junction  is  Qu. 

The  irreversible  thermoelectric  refrigerator  is  composed  of  n 
pairs  of  thermoelectric  elements.  Each  element  is  composed  of 
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P-type  and  N-type  semiconductor  legs.  The  thermoelectric  refrig¬ 
erating  element  is  assumed  to  be  insulated,  both  electrically  and 
thermally,  from  its  surroundings,  except  at  the  junction-reservoir 
contacts.  The  internal  irreversibility  is  caused  by  Joulean  electrical 
resistive  loss  and  heat  conduction  loss  through  the  semiconductor 
between  the  hot  and  cold  junctions.  The  Joulean  loss  generates  an 
internal  heat  I2R,  where  R  is  the  total  internal  electrical  resistance 
of  the  semiconductor  couple  and  I  is  the  electrical  current  generat¬ 
ing  from  the  semiconductor  couple.  The  conduction  heat  loss  is 
K(T'H2  -  ri2),  where  I<  is  the  thermal  conductance  of  the  semicon¬ 
ductor  couple,  T’H2  is  the  hot  junction  temperature,  and  T'n  is  the 
cold  junction  temperature.  Finite  rate  heat  transfers,  i.e.  the  tem¬ 
perature  differences  ( TH2  -  T'H2)  and  (Vl2  -  TL2),  where  TH2  and 
Ttz  are  the  temperatures  of  the  heat  sink  and  heat  source  of  the 
thermoelectric  refrigerator,  respectively,  cause  the  external  irre¬ 
versibility.  For  the  thermoelectric  refrigerator,  the  rate  of  heat 
transfer  at  the  hot  junction  is  Qh2,  and  the  rate  of  heat  transfer  at 
the  cold  junction  is  Qj2. 

Assume  that  the  four  heat  exchangers  among  the  hot  and  cold 
junctions  of  the  thermoelectric  refrigerator,  thermoelectric  gener¬ 
ator  and  their  respective  reservoirs  are  counter-flow,  and  the  heat 
conductance  (product  of  heat  transfer  coefficient  and  heat  transfer 
surface  area)  of  the  heat  exchangers  are  kmFm,  kuFu,  kmFH2  and 
kaFL2,  respectively,  where  km,  kLU  kH2  and  kL2  are  the  heat  transfer 


heat  exchangers  of  the  irreversible  combined  thermoelectric  de¬ 
vice  is  finite  and  F  =  Fm+  Fu  +  FH2  +  F,2  holds. 

Assuming  that  the  heat  transfers  among  the  hot  and  cold  junc¬ 
tions  of  the  thermoelectric  generator  and  the  thermoelectric  refrig¬ 
erator  and  their  respective  reservoirs  obey  Newton’s  law  gives: 

dm  =  kmFm  (Tm  ~T'm)=m  [a/TH ,  +  K(T'm  -  T'u)  -  ^  72r]  (1) 

Qii  =  *uFu (Xu  ~  Tn)  =  m  [otlTu  +  K(T'm  -T'n)+ 1  727?]  (2) 

Qm  =  kH2FH2  (T'H2  -  Tm)  =  n  [«/Th2  -  K(tm  -  T'2)  +  1/2rJ  (3) 

0x2  =  k,2Fu (TL2  -  VL2)  =  n  [otiTu  -  K(T'H2  -  ta)  -  ^727?]  (4) 

where  <x  =  aP-  a.N,  <xP  and  xN  are  the  Seebeck  coefficients  of  the  P- 
and  N-  type  semiconductor  legs  for  each  thermoelectric  power  gen¬ 
eration  and  refrigeration  element. 

3.  Performance  analysis 

Combining  Eq.  (1)  with  Eq.  (2)  gives  the  hot  junction  tempera¬ 
ture  Tm,  and  the  cold  junction  temperature  T'u  of  the  thermoelec¬ 
tric  generator: 


0.5ocRm2/3  -  (0.5 m  +  m2K)RI2  +  maT„,/  -  mK(kmFmTm  +  knFnTn)  -  kmFmknFnTm 
m2a2/2  +  moc(kmFm  -  kuFu)I  -  rnK(kmFm  +  knFn)  -  kmFmknFn 

-0.5a Rm2I3  -  (0.5m  +  m2K)RI2  -  maTmI  -  rnK(kmFmTm  +  k,A Fu  Tu )  -  kmFmkuFuTm 
m2a2F+nta(kHiFHr^dafi^ 


coefficients  of  the  four  heat  exchangers,  respectively,  and  Fm,  Fu, 
Fh 2  and  F/2  are  the  heat  transfer  surface  areas  of  the  four  heat 
exchangers,  respectively. 

The  total  number  (M)  of  thermoelectric  element  pairs  of  the 
irreversible  combined  thermoelectric  device  is  finite  and 
M  =  m  +  n  holds.  The  total  heat  transfer  surface  area  (F)  of  the  four 


Combining  Eq.  (3)  with  Eq.  (4)  gives  the  hot  junction  tempera¬ 
ture  T'h 2  and  the  cold  junction  temperature  T'l2  of  the  thermoelec¬ 
tric  refrigerator. 

r  — 0.5aRn2/3  -  (0.5n+n2K)RI2 -naT^I-nKfk^FmT^  +kL2FI2Tl2)-kH2FmkuFI2Tm 

m~  n2a2I2  +na(kuFL2-kH2FH2)I-nK{kmFH2+kuFI2)-kH2FH2kI2FL2 


(7) 


(8) 


-0.5xRn2!3  -  (0.5 n  +  n2I<)RI2  +  naTH2I  -  nI<(kH2FH2TH2  +  kL2FL2Tt2)  -  kH2FH2kL2FL2TH2 
n2a2I2  +  na (kL2FL2  -  kH2FH2)I  -  nK(kH2FH2  +  ki2FI2)  -  kH2FH2kL2FL2 


Substituting  Eqs.  (5)— (8)  into  Eqs.  (1 )— (4)  yields: 

Qm  =  kmFm{Tm  -  [0.5aRm2/3  -  (0.5m  +  m2K)RI2  +  maTmI 
-mK(kmFHlTm+kuFL1Tn) 

-  kmFmkuFuTm]/[m2a2I2  +  rna(kmFm  -  kuFu'jl 

-  rnK(kmFm  +  kuFu)  -  kmFmknFu]}  (9) 

Qi,  =  fcti Fti { [-0.5aRm2/3  -  (0.5m  +  m2K)RI2  -  rnocTmI 
-mK(kmFmTm  +  kLiFaTa) 

-  kmFmknFuTm}/[rn2a2I2  +  rnu(kmFm  -  knFn)I 

-  rnK(kmFm  +  kuFu)  -  kmFmknFu}  -  Tu}  (10) 

0.H2  =  kH2FH2{TH2  ~  [-0.5 aRn2!3  -  (0.5 n  +  n2K)RI2  -  m.TH2I 

-  nK(kH2FH2TH2  +  kL2FaTL2) 

-  kH2 FH2 kL2 FL2 IH2] / [n2 a2 /2  +  m.(kL2FL2  -  kH2FH2)I 

-  nI((kH2FH2  +  kL2FL2 )  -  kmFmknFu]}  (1 1) 

Q12  =  kt2Ft2{[-0.5aRn2/3  -  (0.5n  +  n2K)RI2  +  mTH2I 

-  nK(kH2FH2TH2  +  kL2FL2TL2) 

-  kH2 FH2 kL2 F,2 IH2]  / [n2 a2 /2  +  ny.(kL2Fn  -  kH2FH2)I 

-  nK(k„2FH2  +  kL2FL2)  -  kH2FH2k,2FL2]  -  TL2}  (1 2) 

The  overall  system  is  a  closed  loop  circuit,  and  the  heat  flow  of 
the  system  is  in  balance,  one  has: 

Q.m+Qi2  =  Q.n+Q.H2  (13) 

Substituting  Eqs.  (9)-(12)  into  Eq.  (13)  and  re-arranging  the  re¬ 
sults  yields  the  equation  that  the  system  stable  working  electrical 
current  should  be  satisfied: 


+3kmFmKRkL2FL2  -  2oc2kmFmTmknFu  +  2(X2kH2FmTH2kmFm 
-  ki2FL2Ti2y2knFu  +  kmFmKRkH2FH2  -  3kuFuKRkH2Fm 
+  2a2 knFuTu km Frn  ~  2kH2FH2TH2a2kuFLt) 

+  m2(-kuFuRkH2FH2kL2FL2  -  kmFmRkH2FH2kL2FL2) 

+  n2(kL2FL2RkuFukHIFHI  -  kH2FH2RkuFukmFm) 

+  mn(2kI2F,2RkuFukHIFm  +  2knFuRkH2FH2kl2FL2 
+  2kH2FH2RkuFu  kmFm  +  2kmFmRkH2FH2kL2FL2)}  (1 7) 

A,  =  ~2m2nK(kI2FL2  +  kH2FH2)(kuFuKR 

+  KkmFmR  +  a2kuFL  iTn  +  a2kmFmTH1) 

-  2m  n2K(kmFm  +  knFu)(I<RkL2FL2  +  KRkH2FH2 
+  v2kH2FH2TH2  +  c/2kL2FL2TL2) 

-  2m2kl2Fl2kH2FH2(KRkIAFu 

+  KRkmFm  +  a2  kn  FnFn  +  a2kmFmTH1) 

-  2n2knF11  km  Fm  (KRkL2FL2  +  KRkH2FH2 
+  a2kH2FH2TH2+u2kL2FL2TL2) 

+  2mn[a2  (kmFm  ku  Fu  kH2FH2Tm 

-  kL2FL2kH2FH2kmFmTa  +  kuF,_,kHIFHlkl2FL2TL, 

-  kmFmkuFakuFnTm  -  ku F(1  km Fm kH2FH2TiA 
+  kL2 FL2 kH2 FH2 ku FuTL2  -  kmFmkL3FnkL2FL2Tm 

+  kH2FH2kL2FL2kmFmTH2)  -  RK(kH2FH2kL1FukmFm 
+  ki2  Fu  kH2  Fh2  kmFm  +  kL2FL2kuFnkHIFHI 
+  knF,  i/<l2F(2kH2FH2)]  -  2mRkHI  Fm  ku Fu  kH2FH2kl2FL2 

-  2 nRkL2FL2kH2FH2kuFnkmFm  (18) 


A4I4+A3I3+A2I2+A1I  +  Ao  =  0 


(14) 


A4  =  oc3m2n2R(kuFu  -  kHI Fm  -  kL2FL2  +  kmFm)  (15) 

A3  =  a2mn[mn(2kH2FH2KR  +  2kH2FH2TH2y2  +  2kL2F,2KR 

+  2kmFmKR  +  2  kL2FL2Tl2a2  +  2a  2k11FL1T11  +  2  knFnKR 


+  kl2Fl2RkuFn  -  k,2Fl2RkmFm  +  2kI2F,2kH2FH2R) 

+  n(-kH2FH2RkuFu  +  2k„1F„1kllF11R  -  kL2FL2RkmFm 
+  k,2F[2RkuFu  +  kH2FH2RkmFm )]  (1 6) 


A2  =  a[m2n(-3kuFuKRkH2FH2  -  2v2kuFuTukH2FH2 
+  2a2kiiFn  T^I<i2Fi2  +  2  y2  kmFm  TH1  k!2FL2 
+  kuFuKRkL2FL2  -  2y2 km FmTm kH2 FH2 

-  2kL2Fl2TL2y2kH2FH2  +  2 k„2FH2TH2a2k12F12 

-  kmFmKRkH2FH2+3kmFmKRkL2FL2) 

+  mn2(—knFnKRki2  FL2  +  2u2kL2FL2TL2kmFm 


Aa  =  2 rnn(KkuFnkmFmTmkH2FH2  -  KkL,FukmFmTukH2FH2 
+  KkLlFukmFmTmkaFa  -  KkuFnkmFmTnkL2FL2 

-  Kk,2Fl2kH2FH2TH2kL]Fu  +  Kkl2F,2kH2FH2Tl2kmFm 
+  KkL2FaTL2kH2FH2kuFu  -  KkL2FL2kH2FH2TmkmFm) 

+  2rn(kmFmknFnkH2FH2kL2FL2Tm 

-  knFukmFHrkmFmkaFaTLi) 

+  2n(kL2FL2TL2kH2FH2kLiFnkmFm 

-  kmFmkaFuhiFukmFmTm)  (19) 

For  the  given  parameters,  one  can  obtain  four  theoretical  solu¬ 
tions  by  solving  Eq.  (14).  Analysis  shown  that  there  is  only  one 
solution  Is,  which  satisfies  /  >  0,  Qhi  >0,  Qu  >  0,  Qh2  >  0,  and 
Q12  >  0.  Upon  that  one  has  the  cooling  load  and  COP  of  the  com¬ 
bined  thermoelectric  device  as  follows: 

Qu  =  ki2F  i2  { [-0.5aRn2/3  -  (0.5n+n2K)RI2 + nyTH2Is  -  nI((kH2FH2TH2 
+  kL2FL2TL2)  -  kH2FH2kL2FL2Tm}/[n2a2I2  +  mx(kL2FL2  -  kH2FH2)Is 

-  (kmFm  +  kL2FL2)  -  kH2FH2kI2Fl2}  -  TL2}  (20) 
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e^Qu/Qm 

kaFn{[-0.5aRn2I3s-(0.5n+n2K)Rl2s+naT„2ls-nK(kmFmTm+kaFaTu)-kmFmkL2FL2TH2\/[n2x2^+na(kI2FL2-kmFH2)h-nK(kH2Fm+knFt2)-kH2FmkL2Fa]  -ru} 
~kmFmljm-[0.5aRm2I2-(0.5m+m2K)KI2+maTmIs-mK(kmFmTm+kaFuTn)-kmFmkuFaTm]/[m2x2I2+ma(kmFm-kaFL,)Is-mK(kmFm+kuFu)-kmFmknFL^} 

(21) 


Obviously,  besides  the  performance  parameters  of  thermoelec¬ 
tric  element  pairs  (a,  R  and  K),  parameters  concerning  the  heat 
transfer  (km,  kIA,  km.  kL2,  Fm,  FM,  F,l2  and  FL2)  would  have  influ¬ 
ences  on  the  device  performance. 

If  kmFm  =  fej.i F(1  =  kH2F,2  =  kL2Fu  ->  oo,  Tm  =  Tm,  Tu  =  T'n, 
Th 2  =  T'h 2  and  TL2  =  T'(2,  Eqs.  (20)  and  (21)  become  the  results  of 
conventional  non-equilibrium  thermodynamic  analysis  [26,27], 


4.  Numerical  examples 

4.1.  Design  variables  and  fixed  parameters 

For  independent  thermoelectric  generator  and  thermoelectric 
refrigerator  devices,  there  is  an  optimum  heat  transfer  area  alloca¬ 
tion  between  the  high  temperature  side  and  the  low  temperature 


Fig.  2.  Working  electrical  current  (a),  cooling  load  (b)  and  COP  (c)  vs.  the  ratio  of 
number  of  thermoelectric  elements.  The  solid  line  represents  result  obtained  by 
using  the  combination  of  finite  time  thermodynamics  and  non-equilibrium 
thermodynamics.  The  dotted  line  represents  result  obtained  by  using  non-equilib¬ 
rium  thermodynamics,  i.e.  effects  of  heat  transfer  are  not  taken  into  account. 


M 


Fig.  3.  Maximum  cooling  load  (a),  maximum  COP  (b)  and  optimum  working 
electrical  currents  (c)  vs.  total  number  of  thermoelectric  elements.  In  (a)  and  (b),  the 
solid  line  represents  result  obtained  by  using  the  combination  of  finite  time 
thermodynamics  and  non-equilibrium  thermodynamics  while  the  dotted  line 
represents  result  obtained  by  using  non-equilibrium  thermodynamics,  i.e.  effects 
of  heat  transfer  are  not  taken  into  account.  In  (c),  the  solid  line  represents  optimum 
working  electrical  current  for  maximum  cooling  load  and  the  dash  dot  line 
represents  optimum  working  electrical  current  for  maximum  COP. 
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Fig.  4.  Contour  maps  of  the  cooling  load  (a)  and  COP  (b)  vs.  ratios  of  heat  transfer 
surface  area.  /,  is  generator  heat  transfer  surface  area  ratio,  and  f2  is  refrigerator 
heat  transfer  surface  area  ratio. 

side  heat  exchangers  for  the  fixed  total  heat  transfer  area  of  the 
thermoelectric  generator  or  thermoelectric  refrigerator.  For  the 
combined  thermoelectric  device,  besides  optimum  heat  transfer 
area  allocation  between  the  high  temperature  side  and  the  low 
temperature  side,  there  is  an  optimum  heat  transfer  area  allocation 
between  the  thermoelectric  generator  and  the  thermoelectric 
refrigerator  for  the  fixed  total  heat  transfer  area  of  the  whole  device. 

In  order  to  describe  the  allocation  of  the  heat  transfer  area, 
three  ratios  of  heat  transfer  surface  area  are  defined:  total  heat 
transfer  surface  area  ratio /=  Fi/F,  where  Fj  =  Fm  +  F;  l,  i.e.  the  ratio 
of  the  heat  transfer  surface  area  of  the  thermoelectric  generator  to 
the  total  heat  transfer  surface  area  of  the  combined  irreversible 
device:  generator  heat  transfer  surface  area  ratio  =  Fm /F^  i.e. 
the  ratio  of  the  heat  transfer  surface  area  of  the  high-temperature 
side  heat  exchangers  of  the  thermoelectric  generator  to  the  total 
heat  transfer  surface  area  of  the  thermoelectric  generator;  and 
refrigerator  heat  transfer  surface  area  ratio  f2  =  FH2IF2,  where 
F2  =  Fh 2  +  FL2,  i.e.  the  ratio  of  the  heat  transfer  surface  area  of  the 
high-temperature  side  heat  exchangers  of  the  thermoelectric 
refrigerator  to  the  total  heat  transfer  surface  area  of  the  thermo¬ 
electric  refrigerator.  Then,  one  has  Fm  =ff,F,  Fu  =J(1  -/i)F,  Fm  = 
(1  -JY2F  and  FL2  =  (1  — /)(1  -f2)F. 

In  order  to  describe  the  allocation  of  the  thermoelectric  element 
pairs,  one  ratio  of  numbers  of  thermoelectric  element  pairs  is 
defined:  x  =  m/M,  i.e.  the  ratio  of  number  of  thermoelectric  element 
pairs  of  the  thermoelectric  generator  to  total  number  of 
thermoelectric  element  pairs  of  the  combined  irreversible  device. 
Then,  one  has  m  =  xM  and  n  =  (1  -  x)M. 


F  (m) 

Fig.  5.  Maximum  cooling  load  (a),  maximum  COP  (b)  and  optimum  working 
electrical  currents  (c)  vs.  total  heat  transfer  surface  area.  The  solid  line  represents 
optimum  working  electrical  current  for  maximum  cooling  load  and  the  dash  dot 
line  represents  optimum  working  electrical  current  for  maximum  COP. 


Obviously,  value  ranges  of  the  variables  x,fff,f2  are  [0,1], 

Numerical  calculations  are  performed  in  order  to  analyze  and 
optimize  the  performance  of  the  thermoelectric  refrigerator  dri¬ 
ven  by  a  thermoelectric  generator.  In  the  calculations,  Tm  = 
450  K,  Tu  =  300  K,  TH2  =  300  K,  TL2  =  290  K,  km  =  60  W/K,  ku  =  15 
W/K,  kH2  =  240  W/K,  kL2  =  120  W/K,  a  =  2.1  x  104  V/K,  K=  1.6  x 
10  2  W/K,  and  R  =  1.2  x  10  3  Cl  are  set  [53], 

4.2.  Effects  of  total  number  and  allocation  of  thermoelectric  elements 

The  system  working  electrical  current  (/),  cooling  load  (Qu)  and 
COP  (e)  vs.  ratio  of  numbers  of  thermoelectric  elements  (x)  are 
shown  in  Fig.  2  by  solid  lines,  respectively.  In  the  calculations, 
M  =  200, F=  1  m2,  /=  0.6,  /,  =  0.7,  and  f2  =  0.4  are  set.  There  exist 
two  optimum  working  electrical  currents  corresponding  to  maxi¬ 
mum  cooling  load  (Qu.max)  and  maximum  COP  (£max).  The  maxi¬ 
mum  cooling  load  and  maximum  COP  vs.  the  total  number  of 
thermoelectric  elements  (M)  are  shown  in  Fig.  3a  and  b  by  solid 
lines,  respectively.  In  order  to  compare  the  results  obtained  by 
using  the  combination  of  finite  time  thermodynamics  and 
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non-equilibrium  thermodynamics  with  these  obtained  by  using 
non-equilibrium  thermodynamics,  the  non-equilibrium  thermody¬ 
namic  results,  i.e.  the  effects  of  heat  transfer  are  not  taken  into  ac¬ 
count,  are  also  shown  in  Figs.  2  and  3  by  dotted  lines.  One  can  see 
that  the  effects  of  heat  transfer  are  obvious  and  should  be  consid¬ 
ered  in  the  performance  analysis  and  optimization  of  the  combined 
irreversible  thermoelectric  devices. 

The  optimum  working  electrical  current  (/opt  Qu )  corresponding 
to  maximum  cooling  load  (Qu.max)  and  the  optimum  working  elec¬ 
trical  current  (/opW)  corresponding  to  maximum  COP  (emax)  vs.  the 
total  number  of  thermoelectric  element  pairs  (M)  are  shown  in 
Fig.  3c. 

One  can  see  that  the  optimum  working  electrical  current  (I)  is 
smaller  and  no  longer  proportional  to  x  when  effects  of  heat  trans¬ 
fer  are  taken  into  account.  I  increases  monotonically  and  the  slope 
of  the  curve  decreases  monotonically  with  the  increase  of  x.  I 
changes  little  when  x  is  near  1 . 

Practical  cooling  load  (Q12)  and  COP  (e)  is  smaller.  The  minimum 
x,  i.e.  the  ratio  of  number  of  thermoelectric  elements  correspond¬ 
ing  the  zero  cooling  load  and  zero  COP  is  bigger  when  the  effects 
of  heat  transfer  are  taken  into  account.  However,  heat  transfer 
has  little  influence  on  optimum  x,  i.e.  the  ratios  of  number  of  ther¬ 
moelectric  elements  corresponding  the  maximum  cooling  load  and 
maximum  COP. 

The  maximum  cooling  load  (Qu.max)  is  smaller  and  no  longer 
proportional  to  M  when  the  effects  of  heat  transfer  are  taken  into 
account.  Qt2,max  increases  monotonically  and  the  slope  of  curve  de¬ 
creases  monotonically  with  the  increase  of  M.  The  maximum  COP 
(£max)  is  no  longer  constant  and  decreases  monotonically  with 
the  increase  of  M. 

Both  the  optimum  working  currents  Iopt  Q_a  and  /„ pW  decrease 
with  the  increase  of  total  number  of  thermoelectric  elements, 
and  lopt.Qa  >  Iopt,e  holds. 

4.3.  Effects  of  total  heat  transfer  area  and  allocation  heat  transfer  area 

Fig.  4  shows  contour  lines  of  the  cooling  load  (Qu)  and  the  COP 
(e)  vs.  the  heat  transfer  surface  area  ratios  /)  and  /2,  respectively, 
with  M  =  100,  F=  1  m2,  and  /=  0.6. 

One  can  see  that  there  are  an  optimum  f  and  an  optimum  f2 
corresponding  maximum  cooling  load  Ql 2, max  or  maximum  COP 
Emax,  respectively.  In  fact,  for  fixed  total  number  of  thermoelectric 
elements  M  and  total  heat  transfer  area  F,  there  are  optimum 
parameters  x,  f  f  and  f2  corresponding  the  maximum  cooling  load 
0i2.max  or  maximum  COP  £max,  respectively. 

Fig.  5  shows  the  maximum  cooling  load  Qt2,max.  the  maximum 
COP  £max,  the  optimum  working  electrical  current  ( hpt.on )  corre¬ 
sponding  to  maximum  cooling  load,  and  the  optimum  working 
electrical  current  (Iopt,£)  corresponding  to  maximum  COP  vs.  the  to¬ 
tal  heat  transfer  surface  area  F,  respectively,  with  M  =  100.  One  can 
see  that  both  the  maximum  cooling  load  and  COP  increases  mono¬ 
tonically  with  the  increase  of  F.  The  optimum  working  currents 
lopt.Qu  increases  with  the  increase  of  total  heat  transfer  area  F  while 
I  opt.,-,  changes  little  and  /„pt,Qu  >  IoPt,e  holds. 

Calculations  show  that  when  total  heat  transfer  area  changes, 
the  optimum  variables  x,  f  and  f2  remain  constant  approx¬ 
imately. 

5.  Conclusion 

A  model  of  internal  and  external  irreversible  thermoelectric 
generator-driven  thermoelectric  refrigerator  is  presented  in  this 
paper  by  using  a  combination  of  finite  time  thermodynamics  and 
non-equilibrium  thermodynamics.  Two  analytical  formulae 
describing  the  cooling  load  vs.  working  electrical  current,  and  the 


coefficient  of  performance  (COP)  vs.  working  electrical  current 
are  derived.  The  performance  optimization  of  the  combined  device 
is  performed  by  searching  for  the  optimum  allocation  of  heat 
transfer  surface  area  of  the  four  heat  exchangers  and  the  optimum 
allocation  of  the  number  of  thermoelectric  element  pairs  based  on 
the  optimization  of  working  electrical  current.  All  the  parameters 
should  be  considered  in  the  design  and  application  of  practice  ther¬ 
moelectric  devices  in  order  to  obtain  the  maximum  economy 
benefit. 

The  results  show  that  when  effects  of  heat  transfer  are  taken 
into  account,  the  working  electrical  current  is  smaller  than  that 
by  non-equilibrium  thermodynamics.  The  cooling  load  is  smaller 
and  no  longer  proportional  to  total  number  of  thermoelectric  ele¬ 
ments  while  the  COP  is  no  longer  constant  and  decreases  monoton¬ 
ically  with  the  increase  of  total  number  of  thermoelectric 
elements. 

For  the  fixed  total  number  of  thermoelectric  elements  and  total 
heat  transfer  surface  area,  there  are  optimum  variables,  i.e.  ratio  of 
numbers  of  thermoelectric  element  pairs  x,  total  heat  transfer  sur¬ 
face  area  ratio  f  generator  heat  transfer  surface  area  ratio  f  and 
refrigerator  heat  transfer  surface  area  ratio  f2,  corresponding  the 
maximum  cooling  load  or  maximum  COP,  respectively. 

Heat  exchangers  are  necessary  in  thermoelectric  generator-dri¬ 
ven  thermoelectric  refrigerator  devices.  There  is  always  external 
heat  transfer  loss  in  the  devices.  Conventional  non-equilibrium 
thermodynamic  analysis  cannot  take  into  account  the  external 
heat  transfer  factor.  The  combination  of  finite  time  thermodynam¬ 
ics  and  non-equilibrium  thermodynamic  gives  a  comprehensive 
analysis  and  optimization  of  the  devices  by  taking  account  of 
external  and  internal  irreversibilities  and  obtains  optimal  perfor¬ 
mance  characteristics  closer  to  practice.  The  results  obtained  here¬ 
in  may  provide  guidelines  for  the  design  and  application  of 
practical  combined  thermoelectric  devices.  Some  experiments  will 
be  made  to  confirm  the  performance  of  the  combined  device  in  the 
later  work  by  the  way  of  Ref.  [53]. 
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